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We present a new method of calculating the rate of neutrino-electron interactions in a strong 
magnetic field based on finite temperature field theory. Using this method, in which the effect of 
the magnetic field on the electron states is taken into account exactly, we calculate the rates of 
all of the lowest order neutrino-electron interactions in a plasma. As an example of the use of this 
technique, we explicitly calculate the rate at which neutrinos and antineutrinos annihilate in a highly 
magnetized plasma, and compare that to the rate in an unmagnetized plasma. The most important 
channel for energy deposition is the gyromagnetic absorption of a neutrino-antineutrino pair on an 
electron or positron in the plasma {yi>e <-> e ). Our results show that the rate of annihilation 
increases with the magnetic field strength once it reaches a certain critical value, which is dependent 
on the incident neutrino energies and the ambient temperature of the plasma. It is also shown that 
the annihilation rates are strongly dependent on the angle between the incident particles and the 
direction of the magnetic field. If sufficiently strong fields exist in the regions surrounding the core 
of a type II supernovae or in the central engines of gamma ray bursts, these processes will lead 
to more efficient plasma heating mechanism than in an unmagnetized medium, and moreover, one 
which is intrinsically anisotropic. 

PACS numbers: ll.10.Wx, 13.10.+q, 97.10.Ld, 97.60.Bw 



I. INTRODUCTION 

Neutrino heating and cooling plays an important role in a variety of astrophysical objects. In core collapse super- 
novae (SNe) neutrinos produced deep within the core of the forming protoneutron star (PNS) are thought to deposit 
some fraction of their energy in a semitransparent region above the surface of the PNS, leading to a neutrino driven 
wind |l| and a robust supernova explosion Although it is currently thought that neutrino-nucleon scattering 
provides the bulk of energy transfer in this region, a significant fraction of the energy and momentum exchange be- 
tween the neutrinos and the matter occurs through neutrino-antineutrino annihilation to electron-positron pairs and 
through neutrino-electron scattering. 

More recently, neutrino-electron interactions have been proposed as an energy deposition mechanism for the central 
engines of gamma-ray bursts (GRBs). A large fraction of the binding energy released during the formation of a 
compact object is emitted in the form of neutrinos and antineutrinos. Both neutron star-neutron star mergers and 
the new "collapsar" models for the formation of GRB fireballs rely on neutrino-antineutrino annihilation to electron- 
positron pairs as a mechanism for the transport of energy from the dense and hot regions to regions of low baryon 
loading where a fireball can form and expand ||. Unfortunately, even for the extremely high neutrino densities 
expected in these systems, the cross section of neutrino-antineutrino annihilation in the absence of a strong magnetic 
field is still quite low, which means that the overall efficiency of conversion of gravitational potential energy to fireball 
energy is also low. Thus, it is very difficult to explain the most energetic of GRBs using neutron star-neutron star 
merger models, and even the collapsar models have some difficulty in depositing sufficient energy to drive the fireball 
through the overlying mantle of the star [Q . 

The astrophysical arguments for the existence of supercritical fields in nature have grown stronger recently with 
the discovery of magnetars ||. In these slowly rotating X-ray emitters, field strengths of up to 4 x 10 15 G have been 
inferred. One may speculate that even stronger fields may be present when such objects are born. This discovery 
negates previous theoretical prejudice against the existence of such strong fields, and suggest further examination of 
the effect that strong magnetic fields may have on neutrino-electron processes. 

In the absence of a strong magnetic field, there are only two allowed types of interaction between neutrinos and 
electrons - neutrino-electron scattering (ve^ <-> ve , ve ± <-> ve ), and electron-positron pair creation and annihi- 
lation (yv <-> e + e~). In the presence of a quantizing magnetic field, electrons occupy definite states of momentum 
perpendicular to the field (Landau levels) and conservation of perpendicular momentum between interacting particles 
is no longer required, as some of the momentum may be absorbed by the field. This allows a number of exotic 
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reactions to proceed, where electrons or positrons jump between different Landau levels while interacting with an 
external neutrino current. These processes are very similar to interactions between photons and electrons in a strong 
field, such as single photon pair creation, which is also forbidden in the absence of a strong field. The additional 
neutrino processes allowed in the presence of a strong magnetic field are absorption of a ncutrino-antineutrino pair 
by an electron or positron {yve^ «-> e^), and the absorption or emission of an electron positron pair by a neutrino or 
antincutrino (ve + e~ «-* v, &e + e~ «-> &). Each of these may have an important effect on the neutrino energy exchange 
opacities in a strongly magnetized plasma. 

A number of authors have shown that the rate of neutrino-electron processes are increased due to the presence of 
a supercritical magnetic field (B > B cr — 4.4 x 10 13 G). Kuznetsov and Mikheev |6j considered neutrino-electron 
scattering and the emission of an electron positron pair by a neutrino propagating through a magnetic field in the 
limit that the electrons and positrons are in the lowest Landau level. They concluded that with a very strong field, the 
neutrino could lose only a small fraction of its energy through this process. While we do not analyse this process in 
detail numerically, our results for the neutrino-antincutrino annihilation processes suggest that the restriction that the 
electron and positron be in the lowest Landau level is too strong. A more general calculation may be performed using 
the results presented in this paper. Benesh and Horowitz j?) calculated the rate of neutrino-antineutrino annihilation 
to electron-positron pairs in a strong magnetic field. They examined two restrictive cases - nearly parallel collisions 
and head-on collisions - and found no significant enhancement of the cross section for these types of interactions. Our 
more general calculation shows that there are regions of phase space with a significant increase in the annihilation 
cross-section, though for moderate field strengths, these regions are small. More importantly, we show that the 
new channels for annihilation which appears in a strong magnetic field - gyromagnetic absorption on electrons and 
positrons in the plasma - are more important over an astrophysically interesting range of phase space. 

Bezchastnov and Haensel || calculated the rate of neutrino-electron scattering in a strong magnetic field by evalu- 
ating the matrix element for the Feynman diagrams shown in figure 1, using the exact wave- functions of the electrons 
in the strong field as external states for the calculation. Considering energy ranges and magnetic field strengths 
appropriate for SNe calculations, they showed that while the overall cross-section for this scattering rate was not a 
strong function of magnetic field strength up to 10 15 G, the strong field introduced significant anisotropies in the rate, 
and could lead to interesting parity violating effects. 

Here we adopt a new approach to calculating the rates of neutrino-electron processes in a strong magnetic field. 
Rather than adopting the direct approach of calculating the cross-sections from the diagrams shown in figure 1, such 
as adopted in |jj and ||, we use a technique of finite temperature field theory (FTFT) which relates the interaction 
rates to the imaginary part of the diagram shown in figure 2 (as the energy of our interactions are low compared to the 
M^-boson mass, we use the Fermi theory to describe our interactions). This technique allows all the neutrino-electron 
interaction rates to be calculated on the same footing, and includes the effects of Pauli blocking. We follow the 
approach of Gale and Kapusta j^] which was developed in the context of calculating the rate of dilepton production 
from hadron interactions in heavy ion collisions. It is based on the fact that the rate of dilepton production is related 
to the imaginary part of the polarization tensor. In general, production as well as decay rates of particles in matter 
can be extracted from the imaginary part of self energies computed by using FTFT fic| l . In our case, the interacting 
particles are neutrinos, not hadrons, and our interactions are electroweak interactions and not due to nuclear forces. 
This approach leads to two important simplifications, in that one obtains the total net rate of a given process, taking 
into account both forward and back reactions. Also, numerical evaluation of the rates is greatly simplified, as one can 
calculate all of the rates virtually simultaneously. 

The plan of this paper is as follows. In the first section, the expression relating the rate of a given neutrino-electron 
process to the imaginary part of the polarization tensor, based on the work of Weldon and Gale and Kapusta 
, is presented. This expression is then used to calculate the rate at which neutrinos and antineutrinos annihilate 
to electron-positron pairs in an unmagnetized vacuum, demonstrating that one obtains the same results as the direct 
calculation. 

In the following section, a strong magnetic field is introduced, and the polarization tensor in a strongly magnetized 
plasma is presented. The imaginary part of the polarization tensor is calculated, and expressions for the rates of 
each of the neutrino-electron processes are given. A detailed numerical evaluation of the annihilation rate is then 
performed, and the rate is shown to reproduce the unmagnetized rate in the limit of a weak magnetic field. 

Finally, the behaviour of the annihilation rate with increasing magnetic field strength is examined, and the impli- 
cations for systems in which significant neutrino heating occurs, such as core collapse SNe and GRBs are discussed. 



2 




V e 
FIG. 1. Feynman diagram for electron-positron pair creation through neutrino-antineutrino annihilation. 
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FIG. 2. This is the diagrammatic representation of the polarization tensor used to calculate the rate of electron-positron 
pair production through neutrino-antineutrino annihilation. The imaginary part of this diagram is related to the rate of the 
process shown in figure 1. 

II. BASIC FORMALISM 

Ultimately, we wish to calculate the rate at which neutrinos and antineutrinos interact with electrons and positrons 
in a strong magnetic held in matter through the processes shown in figure 1. Rather than calculate these rates directly 
through the use of the exact electron wave functions in the magnetic field, such as performed by Bezchastnov and 
Haensel ^] in their evaluation of the neutrino-electron scattering rate, we adopt an alternative approach suggested by 
the use of finite temperature field theory (FTFT). This theory has been pioneered by Weldon |n| and also applied 
by Gale and Kapusta || to calculate the rate at which electron-positron pairs are emitted from the fireball in heavy 
ion collisions. 

According to Weldon the decay and production rate of a boson with energy E at temperature T can be related 
to the imaginary part of the boson self energy II by 

r(£) = -iim{n(£)}, (l) 

where r = Td — T; is the difference of the forward rate of the boson decay Td and the rate of the inverse process Tj. 
The decay and production rates are related by the principle of detailed balance 



r, 



- P E/T 



e^, (2) 



from which we get 



T d (E) = -- i _ e _ E/T Im{IT(£)}. (3) 



Instead of considering the decay of a boson we will investigate first the case of neutrino-antineutrino annihilation to 
an electron-positron pair, vv — > e~e + , in an electron-positron plasma at temperature T e and chemical potential /i e . 
The corresponding annihilation rate is given by equation (||) , where the boson is replaced by the neutrino-antineutrino 
pair. To lowest order the self energy is now given by the diagram of figure 2 which contains an electron loop. The 
four momenta of the neutrino and antineutrino are denoted by Q\ and Q2 and of the electron and positron by P\ and 
P2, respectively. Here we use the notation Qi = (qi, qi) with qi = \m\ and Pi — (Ei, pi) with Ei — y'pf 



to is the electron mass. Then in equation E = q\ + qi holds and the self energy is given by a Lorentz tensor II 



in* • 



3 



Next we want to relate the annihilation rate to the differential rate dR u p/d 3 qid 3 q2 describing the absorption of 
a neutrino-antineutrino pair with momenta qi and q2 in the electron-positron plasma. Usually, this rate would be 
calculated from the matrix element M. of figure 1 according to 

dR vV 1 f d 3 p! f d 3 p 2 r , 

{2t:)H\Q 1 +Q 2 -P 1 -P 2 )Y j {\M\ 2 ) 1 (4) 

i 

where there is a sum over initial states and the Fermi distribution function for the electron is given by 

e(g - Me )/T e + 1 ( 5 ) 

and the one of the positron by replacing /i e by —fx e in equation (|J). The distribution functions in equation (^) describe 
Pauli blocking of the produced electron and positron in the plasma. 

Following Gale and Kapusta Q the differential rate is obtained from the annihilation rate Td by taking into account 
the neutrino current 

,// '"" ' ' M^Im{n^(g 1 + Q 2 )}. (6) 



d 3 qi d 3 q 2 2 qi 2q 2 1 - e- E / T ' 
Here the leptonic tensor for the neutrinos is given by 

M"" = v{q 2 )r(l - 75M<ZlM<Zi) 7 "(l - ls)v(q 2 ) 

sis 2 

= KQ2Q1 + QiQl - {Q1Q2W + iQx a Qw^ v \ (7) 

where we have summed over the neutrino spins. This formula differs from the one for the dilepton production rate 
given in the appendix of ref. Q by the absence of a photon propagator and a different structure of the leptonic tensor 
for massless neutrinos instead of electrons leading to the appearance of the antisymmetric tensor s"^ 1 ' in equation 
(^). Furthermore the two expressions differ by a factor exp(E/T e ) according to (|^) because we are considering a 
decay instead of a production rate. 

The total rate for the production of electron-positron pairs from neutrino-antineutrino annihilation follows by 
integrating the differential forward rate, equation (Q), over the neutrino distribution functions j v ,v{qi) 

*» -fgfrj |^[M*>Ata> - (i - /.(*»(! - (8) 

Here we have included the back reaction e~e + — > vv, which is described by the second term in the square brackets. 
It should be noted that our derivation of the rates at no time assumed that the neutrinos are in thermal equilibrium 
with the electron-positron plasma. For example, neutrinos that escape from the core of a supernova explosion may 
have a roughly thermal distribution with a temperature of T v ~ 5 McV. After leaving the core these neutrinos will 
interact with the surrounding electron-positron plasma which has a typical temperature of about T e ~ 1 MeV. In this 
case the neutrino distribution functions in equation (g) are given by the equilibrium Fermi distributions, equation (js|), 
where the temperature and the chemical potential are replaced by the corresponding values for the neutrinos escaping 
from the core, T v and [i v . 

One may also determine the rate at which a neutrino scatters in a thermal plasma in a similar manner. The 
differential rate at which a neutrino (or antineutrino) scatters from an initial state Q\ to a final state Q 2 , is given 
simply by replacing Qi + Q 2 by Qi — Q 2 in Imjn^,,} and E = q% + q 2 by qi — q 2 in equation (^). The total rate of 
scattering, R ve ± , is then given by 



with 



^r^T^^ M ""^ n '- M >- Q ^ < 10 ' 
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where E = qi — q2- Note that both the electron and positron scattering rates are included within the imaginary part of 
the polarisation tensor in equation (|io| ) . The strength of this method lies in the fact that the correct calculation of the 
imaginary part of the polarisation tensor leads the differential annihilation and scattering rates for all of the separate 
channels that arc allowed within the plasma. Thus, if one includes the effect of the magnetic field in the polarisation 
tensor properly, all of the exotic processes which occur in a magnetic field are included in the rate calculations by 
default. 



A. The unmagnetized rate 

To demonstrate that the two approaches, based either on the matrix element as in equation or on the FTFT 
as in equation (^) are equivalent, we use equation (|g) to calculate the rate of neutrino-antineutrino annihilation to 
electron positron pairs in an unmagnetized vacuum at zero temperature, and show that it reproduces the results of 
the calculation of the process using the matrix element method. 

In the one-loop approximation the polarisation tensor in equation (^J) is given by 

U ^ {K) = „^i T r J |^7 M (cv - c Al5 )G(Ph»(c v - c Al5 )G(P - K), (11) 
where Gp is the Fermi constant, G(P) is an electron propagator, and 

cv =±- + 2sm 2 9 w , ^ = ±1 ( 12 ) 

with 9\y the Weinberg angle, and the plus sign corresponding to electron neutrinos, and the minus sign to muon and 
tau neutrinos. This tensor has been formulated in the limit that the neutrino energies are small compared to the 
W mass, which will always be justified in practice. This allows the use of effective Fermi 4-vertices in the processes 
shown in figure 1 and 2. The polarization tensor in equation (^) may be decomposed into vector- vector, axial-vector, 
and axial-axial parts through 



IP W = c 2 v a^ + IcyCAdf + c\a»l (13) 



with 



•^Tr J ^JL^ G (P)rG(P K), (14) 



< = ^Tr J-^l {^ 5 G(Ph"G(P -K) + 7 "G(P) 7 ^ 5 G(P - K)] , (15) 

and 

< = -^ Tl 7 ^l^GiPh-^GiP-K). (16) 

In this section we will neglect matter contributions to the rate, i.e., we set T e = in equation (^), set the distribution 
function to zero and use only the vacuum part of the polarization tensor, equation (p"l|). This corresponds to neglecting 
the effects of Pauli blocking on the final state particles in equation ( pj). Inserting only the vacuum part of the electron 
propagators, the vacuum part of the polarization tensor, equation (jll|), may be written 

n /»Yjn = -^T f r(cv-c Al5 W + m)r(cv-c Al5 )(P-fi + m) 

[ ' 2 J (2tt) 4 [P 2 - m 2 + iO] [(P - Kf - m 2 + iO] ' 1 ' 

Note that the real part of this tensor must be renormalized due to the singularities in the denominator of the integrand, 
but the imaginary part is finite. 

The imaginary part of IF^ (actually, the anti-hermitian part) is related to its discontinuity across a branch cut. 
The branch cut is related to poles that appear above the pair creation threshold, K 2 — 4m 2 , and it lies along the 
positive K 2 axis from 4m 2 to infinity. One has 
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I m {Il^(K)} = -^I^ (18) 



32tt 2 
with 

= J d 4 P Tr W(c v - c Alh ){f + m)Y{cv - c Alb ){? -$T + m)] 

xS(P 2 — m 2 ) S((P — K) 2 — m 2 ). (19) 
Inserting a delta function into equation ( fl9| ) and relabelling P — ¥ P\, leads to 

1^ = -[ d i P 1 J d 4 P 2 Tr [^(cv - CA7s)(fi + m)^{c v - c A ^){f 2 - m)] 

x5{Pf - m 2 ) <5(P 2 2 - m 2 )<5 4 (P 2 + P x - K). (20) 

The trace in the integrand of equation (EG) may be recognized as proportional to the electron leptonic tensor, L^ v , 
and the rate of pair production, equation (p|), is then given by 

#A =—jW?J (2tt)3 2^ 2i? 2 2(?1 2<Z2 (2 " } '(ft+ft-Oi-*). ( 21 ) 

where 



Jd*P l 5(P 2 - m 2 )= (22) 

has been used. 

Equation ( pT| ) is identical to the rate, equation (4), directly calculated for the process shown in figure 1 by electroweak 
theory, using L 11V M^ V — X^d-^l 2 )- Thus we have shown that the procedure outlined above is identical to the standard 
method of calculating rates. 

For later reference, the rate at which the annihilation process proceeds in the absence of a magnetic field at zero 
temperature and in the limit of a vanishing electron mass, m = 0, is given by 

dR _ 2G M4+4i (Z 2 (? 2 (1 _ cose) 2 j (23) 



d 3 qid 3 q2 37r 
where is the angle between the neutrino and antineutrino. 

III. STRONG MAGNETIC FIELDS 

We may now make a significant generalization of our theory, which justifies its introduction. The essential point is 
that the polarization tensor used in equation (^|) contains all the information about the electronic part of the process, 
without reference to the presence or lack of a medium, or the presence or lack of a magnetic field. All that is required 
to calculate the correct rate is the appropriate polarization tensor for the given plasma. To determine the rate of 
neutrino-electron processes in a strongly magnetized plasma, one must use the polarization tensor for a strongly 
magnetized plasma. Thankfully, for our purposes, the polarization tensors in a strong magnetic field are well known, 
the vector- vector part was calculated in the form used here by Melrose and coworkers [0 , the vector-axial part was 
calculated by Kennett and Melrose |T^| , and the axial-axial part may be calculated in a manner directly analogous to 
both. 

The polarization tensor is calculated in [pd| and |T^| using a vertex formulation of the coupling between an electron 
in a magnetic field and an external 4 current. The various parts of the polarization tensor can be expressed in terms of 
the vector vertex function [r^(A')] M , and the axial vertex function, s[r^(A')] A1 , where q' and q denote the quantum 
numbers of the particles on either side of the interaction (py, momentum, a, spin, n, Landau orbital), and e and e' 
denote the nature of the particles, with a plus sign for electrons (the particles) and a minus sign for the positrons (the 
antiparticles). The various components of the polarization tensor may be written, 

n',n— e'e— ±1 



2ir w- ee a + e'e a > + iO 1 ' 



G 



GleB 



47T ^ ^ 

n',n— e',e— ±1 



dp,, 



{|(e'-e)+en*-e'r^} 



X ' 2^ w -ee + e'e ,+iO 5T ^' (25) 



<{k)=-^ e e 

n ; ,n— ±1 



with 



dp|| {^-e) + ^-e^,} 



^r= E [q'4(k)]"[r^(k)]-, (27) 

CT',<7=±1 

= E 5 { 5 [r^(k)]^[r^(k)]- + tr^(k)]" B [r&(k)r } , (28) 



55^= E 5 [r^(k)]^ 5 [r^(k)]-, (29) 

er' : <T=±l 

and where -fC = (w,k), are the electron distribution functions, and e q denotes the energy of the particle in the 
magnetic field, 

e q = e(p,| , n) = (m 2 + p\ + 2neB) . (30) 

Implicit throughout equations (|24|) to ( p9| ) is the relation e'pj, = ep,, — fc„ . Note also that we are using natural units 
with all physical quantities scaled against the electron mass. 

The polarization tensors given in equations (|24|) to ( p6| ) contain the contributions from both the vacuum polarization 
and the electrons and positrons in the plasma. The vacuum part is given by the term i(e' — e) in the numerator 
of the integrand. The real part of this term is divergent and must be renormalized. We are concerned here only 
with the imaginary part, which is finite. Hence, we perform no renormalization here. The electrons and positrons of 
the plasma contribute to the polarization tensor through the distribution functions (more properly, the occupation 
numbers) n q , which in a magnetized thermal medium are given by the Fermi distributions 

ni = n e (p\\,ri) = — r - t . 1 , /rr . (31) 

The general form of the magnetic vertex functions are given in |l2| , and a general discussion of the properties of 
the vector part may be found in p"l| ] . The components of the vertex functions are 

[r^(k)]° - C g ,C q [8 a , i(T (5 e ,, e 13 + aS^ e A)(4_i + Pn'PnJ 1 ^) 

-eS^-riaS^tA - 5 e ,,- e B)(-p n J l +_°_ <7 + p n , j\, _,_„)], (32) 

[r^(k)] 1 = C q> C q %>, a e{5 e ,. £ V + a8 e , ,_ e C)(- p n j\?-l-<, - P»'4-l-*,) 

-<^,_ ff (a<5 e /, e C - *ei,_ e 23)(Jif,_, - Pn'PnJ\^ 2a )l (33) 

[r^(k)] 2 = iC q> C q [5 a ,, a e{o5 fJ , t V + 5^. t C){p n Jl+^_ a ~ p n >4-i + „) 

-6 a ,^ a (6 e ,, e C - <7^,-«X>)(4-J + Pn'PnJfa-to)], (34) 
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and 



[r^(k)] 3 = CtfC q [8 a , ia (6 t ,, e A + a6 e ,,- e B)(4_ l +p n ,p n Jl%) 

-eS^^aS^B ~ 8 t , ^A){- Pn j\t^ a + Pw4-i-.% 

5 [r^(k)]° = C ql C q [8 a ,, a {a8^ e A + ^,_ e B)(4_ ; - Pn'PnA^-i) 

-e8 a <- cr (-8 e <,eB + cr8 < ,i,-eA){-p n Jp~° l _ a - p n '4-i-<r) 



(35) 



(36) 



5 [r^(k)] 1 = C q> C q [6 a ,^(a5^ e C + 5 t ,^V)(- Pn J l +_V a + Pn'4-i+a) 



(37) 



5 [r^ g (k)] 2 = iC q ,C q [8*',cr<8 el , e C + a8 e ,,- e V)( Pn J l +_°_„ + Pn'4-l+a) 
-8 a ,^ a {-o8 t *. t V + 8 e .^C){4_ l - Pn'PnJ^^i 



(38) 



5 [r^(k)] 3 = C q ,C q [8 a ,, a (a8^£ + S el ,- e A){4_ t - PwPnJ 1 ^) 

-8*,-<,e(-8s,eA + a5 t> ,- t B){-p n J l +° l _ (J - p*4-i- a ) 



where I = n — i(l + er) is an orbital quantum number, and 



C q ^ 



\s q +s^(s q+m ) \ 1/2 
le q e q 



(39) 



(40) 



P\\ = P\\/( £ q+ e q), p n =p n /{e° q +m) 7 



(41) 



A = p\+ Ph B = l + p\ P \\ 



(42) 



C = p'\\-Pb V = l -P\\P\\ 



(43) 



and p n — \/2nB, = e(0,n), and the primed versions of the same quantities have unprimed variables replaced by 
primed variables. Note again that e'p'^ — ep\\ — fc|| is implicit throughout. 

The functions J™(kj_/2B) in equations ( |32| ) to ( p9| ) (the argument is suppressed) are generalized functions related 
to the Laguerre polynomials and defined by 



j:( X ) 



1/2 



(n + v)\_ 

where L^(x) are associated Laguerre polynomials. One may write explicitly 

J«(x) = exp(- a; /2)x«/ 2 £ + {n w V)l]1/ l 
* — ' (v + m) (n — mm 

Note that J™ (a;) = if n < or v < —n. A summary of the properties of the J functions may be found in [fl~Tf] . 



(44) 



(45) 



A. Annihilation rate 

We now use the imaginary parts of equations (|24|) to ( [26] ) to determine the rate at which a neutrino annihilates on 
an antineutrino in a strongly magnetized electron-positron plasma. The imaginary parts of these integrals come from 
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the poles in the denominators of the integrands, where lu — ee q + e'e q i = 0. The contribution of these poles may be 
calculated through the Plemelj formula 



= V- 



LU — LUq 



LU — LUq 



— iir5(io — wq), 



(46) 



where V denotes the principle value. The integral over the delta-function introduced into equations (|24|) to (|2q ) by 
the Plemelj formula may be done explicitly by determining the positions of the poles of the integrand. For lu > 0, 
which is the case for neutrino-antineutrino annihilation, the denominators of equations (|24|) to (|26| ) only have zeros if 
e = 1, e' = — 1, (pair creation), e = —1, e' = —1 (neutrino pair absorption on a positron), and e = 1, e' = 1 (neutrino 
pair absorption on an electron). Each of these processes occurs in a different region of the phase space, determined by 



the existence of the roots of the energy conservation equation lo — ee q + e'e q > — 0. Writing S 
has general solutions of the form 



fcij, this equation 



P\\ 



ith 



2 ^ 2 



2(n-n')B 



(47) 



(48) 



and 



A = e 



4(1 + 2nB) 



(49) 



The roots given by equation ([47]) are real if A > 0, which is true if and only if 

2 



< s < 



Vl + 2n'B - V1 + 2nB 



which corresponds to a pair absorption on an electron or positron, and 

2 



<5> 



Vl + 2n'B + Vl + 2nB 



(50) 



(51) 



which corresponds to an electron positron pair creation. Note that for two annihilating neutrinos one has S — lu 2 — k?, > 
0. 

The contributions of the vector-vector, axial-vector and axial-axial parts of the rates may be separated out through 



dR v 



1 



1 



d 3 qi d 3 q 2 2 qi 2q 2 1 - e - E / T - 



M^lm{c 2 v a^(K) + 2c v c A af {K) + c\ a$(K)} 



(52) 



with K = Qi + Q 2 - Each of the polarization tensor components in equation ( |52| ) contains a sum over two sets of 
Landau levels evaluated at the root of the denominator, and a sum over the e and e, which denote the type of particles 
which are interacting. Here, we divide the rate into the three contributing combinations of e and e', corresponding to 
three different physical processes. 



1. Neutrino pair creation 

If equation ( |5l| ) holds, and e = 1, e' = — 1, the poles in the integrand of the polarization tensor components 
correspond to the physical process of neutrino-antineutrino annihilation to form an electron-positron pair, vv <-> e + e~ . 
Letting X stand for a blank, 5, or 55, the imaginary part of the polarization tensor components may be written 



GpCB ^ £q£q< 



x T^(-l + n e -(pu,n) +n e +{p\,,n')) (53) 

Ph :u> — e„ — e„i=0 
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with pj| = — p|| + fc|j. This equation is combined with equation (|32|) and equation (g) to obtain the net rate of 
annihilation to neutrino antineutrino pairs to electron positron pairs. The sum in equation ( p3| ) is evaluated at the 
parallel momenta for which the resonance condition is satisfied 

P|| = ^±|a 1/2 , (54) 

pj| = ^T|A 1/2 . (55) 
Equation should be evaluated at both of these roots, and the sum over the Landau orbitals is limited to n and 



n' such that y/l + 2nB + VI + 2n'B < ^lo 2 - fcjj. 

The form taken by the blocking factors in equation ( |53|) is due to the fact that the imaginary part of the polarization 
tensor is related to the difference of the forward and reverse rates jl(| . 

2. Gyromagnetic absorption on positrons 

If condition ( |50| ) holds and if n' > n, with e = — 1, and e' = — 1, the neutrino- antineutrino pair may be absorbed 
on a positron, vve + <-> e + , a process which is forbidden in the absence of a strong magnetic field. The contribution 
to the imaginary part of the polarization tensor due to the root corresponding to this process may be written 

lm{c%(K)} = ^?- Y -. £q£q ' , , x 

xT^_(-n e +(p»,n) +ra e +(p(i,n')) (56) 

PH-w+e q -e q /=0 

with p'| = p|| + fc|| . The sum is evaluated for this process at parallel momenta 



h 

2 * 2' 



Pn = -4e±^ 1/a , (57) 



pj,=^±^A 1/2 . (58) 



3. Gyromagnetic absorption on electrons 

On the other hand, if condition ( pOj ) holds and if n > n', and e = 1 and e' = 1, the neutrino-antineutrino pair may 
be absorbed on an electron, vve~ <-> e~, a process which is also forbidden in the absence of a strong magnetic field. 
The contribution to the imaginary part of the polarization tensor due to the root corresponding to this process may 
be written 



ri2 r> 00 



xT^_(n e - (p\\ , n) - n e - (p'« , n')) (59) 

p\\:u-s q +e q i=0 



with p'| = p|| — fc|| . The sum in equation (|59| ) should be evaluated at the roots 

fen . u 



P[| = ^±^A 1/2 , (60) 



P\l=-^±U 1/2 - (61) 
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For the gyromagnetic absorption processes, the limits on the sums over the Landau levels enumerated by n and 

n' are given by \y/i + 2nB — \/l + 2n'B\ > y ^ lo 2 — fcjj. There is no upper limit on either n or n' implied by this 

relation. Instead, the upper limit on the Landau level is constrained by the fact that one of the distribution functions 
in equations ( p6| ) and ( p9| ) must be non-zero for there to be a rate, and these distribution functions are sampled at 
the resonant momenta given above, which grow with increasing Landau number. 



B. Scattering rate 

A similar procedure may be used to calculate the rate of scattering of a neutrino or antineutrino off a thermal pair 
plasma in a strong magnetic field. In this case, the frequency of the disturbance in the field is given by lo = q\ — q 2 , 
which may be negative or positive. Hence, there are four contributions to the scattering rate, whereas there are only 
three contributions to the annihilation rate. The processes correspond to e = l,e' = — 1, pair creation by a neutrino, 
e = l,e' = — 1 pair absorption by a neutrino, e = — l,e' = —1, neutrino positron scattering, and e = l,e' = 1 neutrino 
electron scattering. The rates of these processes are given by equation (||) with 

^ - 2 gl 2 g2 l-e-^ M ^4^W + 2c v c A aTiK) + c\ (62) 

and with K = Q\ — Qi- Again, the scattering and absorption processes occur in different regions of the phase space, 
with the conditions 



hi 



< 



VI + 2n'B - v 7 ! + 2nB 



(63) 



which corresponds to neutrino scattering off electrons and positrons, and 

2 



5> 



v 7 ! + 2n'B + VI + 2nB 



(64) 



which corresponds to electron positron pair creation and absorption on neutrinos. 



1. Electron-positron pair creation and absorption 



The contribution of the pair creation and absorption processes to the total neutrino scattering rate, ve + e <-> v, 
appears in two parts, for opposite signs of the energy transfer, lo. For to = qi — q 2 > 0, one has a contribution for 
e = 1, e' = -1 ii 8 > [Vl + 2nB + Vl + 2n'B] 2 . This contribution is given by 

G 2 e eB £q£q' 



xT^"_ (-1 + n e - , n) + n e+ , n')) 



p«:u-e g -e ,=0 



with p'| = — p|| + /c|| . Equation ( |65| ) should be evaluated at both the roots given by, 

p„4^ Al/2 ' 



(65) 



(66) 



Pil = ?^~A 1/2 



(67) 



For negative energy transfer, lo = (71 — 92 < 0, one has a contribution for e = —1, e' = 1 if <5 > 
[vT + 2nB + \/l + 2n'B\ . This contribution is given by 



Im{<(^)} = 



GleB 



4tt 



E 



e„e 



qt q > 



X 



Tr + (1 - n e + (p|| , n) - n e - (py , n )) 



P||:a;+e,+e g /=0 



(68) 
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with p'| = — p\\ — k\\ . Equation ( |68| ) should be evaluated at both the roots given by, 



P ||=-^±|A^, (69) 
Pl,=-^^A 1/2 (70) 



2. Neutrino-electron scattering 

There is a contribution to the imaginary part of the polarization tensor for e = e' = 1, if 8 < 
[vT+ 2nB + vT+ 2n'B\ . That is, if this condition holds, the neutrino-electron scattering is kinematically allowed, 
and it contributes to the total neutrino scattering rate through 

^2 



lm{<(*0} = - 1 



GieB 



An n^=M £ l'- p '\\ £c > 



xT++K- {p\\ , n) - n e - , n')) (71) 

p n :uj-e q +e q ,=0 

with p'| = p|| — fcy . If 5 < 0, then there is only one root to the equation uj — e q + s q i = 0, which lies at 

P|| = ^ + |A^, (72) 

if fc|| > 0, and at 

P|| = ^-|AVa (73) 

if fc|| < 0. If, on the other hand, < S < [\ZT~\~ 2nB + Vl + 2n'B\ 2 , then there are possibly two roots or no roots, 
depending on the relative sign of uj + fcy and n' — n. If n' — n has the same sign as ui + k\\, then there are no roots, 
and this scattering reaction is not kinematically allowed. If they are different in sign, then both the solutions given 
by equations (J73) and (|73[) are roots to the equation and represent open scattering channels. 



3. Neutrino-positron scattering 

The structure of the theory for neutrino-positron scattering is similar to that of neutrino-electron scattering. The 
rate of neutrino positron scattering may be written, 



4n n pio\P\\ e «'-P'f«\ 
X T^_ (-n e+ (p|| , n) + n e+ (pj, , n')) (74) 



with p',, = + fc|| . If lo 2 — fcn < then the energy conservation equation, u> + e q — e q > =0 has a single root given by 



k 



P\\ 



-^ + ^, (75) 



2*2 



if fey > 0, and at 



if fen < 0. 



Alternately, if < 6 < U/l + 2nB + y/T+ 2n'S] 2 , then there are possibly two roots or no roots, depending on the 
relative sign of lu + fey and n' — n. In this case, if n' — n has the opposite sign to u> + fen , then there are no roots, 
and this scattering reaction is not kinematically allowed. If they have the same sign, then both the solutions given 
by equations d73) and ([76h are roots to the equation and represent open neutrino-positron scattering channels. 
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IV. EXPLICIT EVALUATION OF THE RATES 



In equations ( |53| ) to (|59j) and equations (|6£]) to (|74|) we have expressions for the rates of the various neutrino- 
electron processes in terms of double sums over the electronic Landau orbitals. In general, these expressions may only 
be evaluated numerically. To demonstrate how this is done, we now make a simple calculation of the rate of electron- 
positron pair production through neutrino- antineutrino pair annihilation in such a strong field that only electrons 
and positrons in the lowest Landau orbital may be produced. This is not a physically interesting calculation, as the 
approximations made to allow an analytic calculation are quite extreme. 

A. Emission into the lowest Landau orbital 



To illustrate the use of equation ([32|) we calculate the rate at which electrons and positrons are created in the lowest 
Landau orbital by neutrino-antineutrino annihilation for head-on collisions. 

We assume that the magnetic field lies in the ^-direction, and let the neutrino 4-momenta be given by Q\ = 
(<7i, gi sin a, 0, gi cos a) which lies in the x-z plane, and Q2 — (92, 1i sin 6* cos 4>, q 2 sin sin <jf>, qi cos 9). As the electron 
and positron are created in the lowest Landau orbital, we have that n = n' = 0, I = I' = 0, and thus, a = a' = — 1. 
This then leads to the following expressions for the vertex functions, 

r° = -C q C' q (p\ l + p ll )Jl (77) 
Y^ = -C q C' q {l + p\ n )J° , (78) 
b T» = C q C' q (l + p\p\\)Jl (79) 
5r 3 = C g G>j |+(O ||)J (80) 

with the other components zero. 

Using equations ([77]) to ( |S0| ) , the contraction over the neutrino tensor in equation ( |52| ) may be performed explicitly, 
leading to the rate of production of electron positron pairs in the lowest Landau orbital from a neutrino-antineutrino 
pair, 

dRvv _ G 2 F eB ^^'2 j^-k 2 t I2B £ q £ q' 
<P qi <Pq 2 2ir q q " \ Pll e q > -p'f q \ 

x [(1 + cos 9 cos a + sin 9 sin a cos <j))(cvA — c^B) 2 
+2(cos# + cosa)(c v A — caB)(caA — CyB) 

+ (1 + cosflcosa — sin 9 sin a cos (f>)(cvB — caA) 2 ], (81) 

where A = p\\ + p',,, B = 1 + p\\p'^, and (Jq) 2 = e ~ fc J-/ 2S has been used. 

Equation ( |Sl| ) is to be evaluated at both of the solutions to u ~ e(j>\\, 0) + e(— p\\ + fey, 0). These are given by 

(82) 




with 

u = qi+q2, fcjj = q\ cos 2 a + q\ cos 2 9 + 2qiq 2 cos a cos 6*. (83) 

Also, 

h.. , , / A 

(84) 
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We further restrict our attention to head-on collisions between a neutrino and antineutrino with the same energy, 
q± = Q2 = q. Leaving the angle between the neutrino and the magnetic field as a, the antineutrino angles must be 
9 = 7r — a, and <f> = ir. Under these conditions, fcy = 0, and k± = 0, as we arc in the centre of momentum frame. 

Thus, p|| = ±-\/g 2 — 1, and p'u = —p\\- Evaluating equation ( |Sl| ) leads to 

dR vP G%cleB m 2 e . 2 

— — — : :Sin a, (85) 



d 3 q 1 d?q 2 tt qy/qi - 1 

where there is a factor of 2 due to the two roots of the resonant denominator. Note the dependence on the angle 
with respect to the magnetic field. This shows the intrinsic anisotropy of this process, in that neutrinos approaching 
head-on along the field line cannot interact to produce pairs. 



B. General Calculation 



In general, the calculation of the rates proceeds numerically. Our simple calculation motivates the numerical 
approach. One specifies the incident neutrino energies and angles, and explicitly calculates the leptonic tensor, 
equation (Q). Then, for each allowed n and n' for the process in question, one calculates the location of the roots, 
and the vertex functions, equations ( |32] ) to (|3^), evaluated at these roots. From the vertex functions, one constructs 
the xTjf^ matrices and contracts them with the leptonic tensor. The results of the sum are then premultiplied with 
the phase space factors and other prefactors to give the differential rate of the process. 

While this is a reasonably complicated procedure, it has a certain natural advantage, as it exploits the similarities 
between the calculations of the different processes. The alternate method of Bezchastnov and Haensel || is not 
essentially simpler. The procedure outlined here can be performed in such a way that the rates for all seven processes 
can be calculated simultaneously, with the common elements of the calculation shared for computational efficiency. 



C. Weak magnetic field limit 



Annihilation 

10 1 I 1 I I 1 I r 




FIG. 3. Rate of neutrino-antineutrino annihilation for a nearly head-on collision between a 5 MeV neutrino and a 6 MeV 
antineutrino in vacuum as a function of magnetic field. The rate has been scaled against the unmagnetized rate, and shows 
that the unmagnetized theory is reproduced in the limit of a small magnetic field. The behaviour of the rate for very high field 
strengths where only particles in the lowest Landau orbitals are produced is linear, and follows the results derived in section 
4A. 

Thus far, it has proved impossible to show analytically that the expressions given above reduce to the unmagnetized 
rate in the limit of a weak magnetic field. However, the numerical calculation of the rates does reproduce the expected 
behaviour in the weak field limit. In figure 3, we show the rate of electron-positron pair production through neutrino 
annihilation as a function of magnetic field strength, normalized to the unmagnetized rate. This plot shows that the 
unmagnetized limit is obtained in the limit of small magnetic field strength, as the graph tends to unity for a small 
magnetic field. 
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D. Numerical results 



We turn now to a more general survey of the properties of the annihilation rate as a function of magnetic field 
strength, neutrino energy, plasma temperature, and the various angles in the system. Throughout we assume that the 
chemical potential of the plasma is zero. In general, the rates are very complicated functions of these variables. To 
determine the effect of a strong magnetic field on an astrophysical system, these rates must be averaged over thermal 
distributions of neutrinos and over a realistic angular distribution. This exercise is beyond the scope of this work. 

In figure 4 we show the rate of neutrino-antincutrino annihilation, scaled to the unmagnctized rate, of a 4 MeV 
neutrino and a 5 MeV antincutrino in a pair plasma of temperature T = 2 MeV as a function of magnetic field 
strength. This plot shows the fine structure that comes about due to the quantized perpendicular momentum states. 
The underlying rate of electron-positron pair production is a weak function of the field strength until around 10 15 G. 
This is the field strength at which B w q\ + in scaled coordinates, and thus, the magnetic field begins to dominate 
above this level. It is also here that the rate of gyromagnetic absorption of neutrino- antineutrino pairs becomes 
significant - this is because the amount of energy required to lift an electron or positron from its lowest Landau level 
to another (low) Landau level is more closely matched to the available neutrino energy. This allows the phase space 
for this reaction to grow and for the process to proceed more vigorously. Consequently, the effect of the magnetic 
field is stronger for low energy neutrinos. 

In figure 5 we show the annihilation rate for the same conditions as figure 4, except that the magnetic field is fixed 
at 4 x 10 15 G and the energy of the annihilating neutrino is allowed to vary. This plot shows the standard increase in 
cross section as a function of energy for the pair creation rate, but also shows a number of other interesting features. 
For instance, for low energy neutrinos, the pair creation rate is suppressed - this is because there is only enough 
energy to create pairs in the lowest Landau level, and these states are filled by the thermal plasma. On the other 
hand, there is a significant rate of gyromagnetic absorption on the plasma particles in the lowest Landau orbital, even 
for annihilation with the lowest energy neutrinos. Of course, these rates must be averaged over thermal distributions 
of neutrinos and antineutrinos, so it is very unclear as to which of the processes will dominate for which regions of 
parameter space. 

Figures 6 to 8 show the annihilation rates as a function of the angles of the system. They show that the annihilation 
rates are anisotropic, and can have a very sensitive dependence on angle. At higher neutrino energies or at lower field 
strengths there can be many more resonances, leading to highly anisotropic heating and momentum deposition. 



Annihilation rate in a 2MeV plasma 




Magnetic field strength [G] 



FIG. 4. Representative calculation calculation of the rate of annihilation of a 4 MeV neutrino and an 5 MeV antineutrino in 
a 2 MeV pair plasma as a function of magnetic field strength. The solid line shows the rate at which the pair process occurs as 
a function of field strength, and the dashed line shows the rate of gyromagnetic absorption onto electrons and positrons. The 
angular variables were chosen such that a = 0.2, 6 = 1.1, and 4> = 0.4. 
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FIG. 5. The rate of neutrino-antineutrino annihilation in a 2 MeV pair plasma as a function of the neutrino energy. The 
magnetic field is fixed at 4 x I0 15 G, and the other variables are as in figure 4. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

alpha angle 

FIG. 6. The rate of neutrino-antineutrino annihilation in a 2 MeV pair plasma as a function of a, the angle between the 
incoming neutrino and the magnetic field. The other variables are as in figures 4 and 5. 




0.0 0.5 1.0 1.5 2.0 2.5 

theta angle 



FIG. 7. The rate of neutrino-antineutrino annihilation in a 2 MeV pair plasma as a function of 6, the poloidal angle between 
the incoming anti-neutrino and the magnetic field. The other variables are as in figures 4 and 5. 



1G 



Annlhilotion rote 



| | | | | , 

12 3 4 5 6 

phi angle 

FIG. 8. The rate of neutrino-antineutrino annihilation in a 2 MeV pair plasma as a function of <f>, the toroidal angle between 
the incoming anti-neutrino and the magnetic field. The other variables are as in figures 4 and 5. 



V. CONCLUSIONS 



We have presented here a new calculation of the rates of neutrino-electron interactions in a strongly magnetized 
thermal electron-positron plasma using FTFT. Starting from the imaginary part of the polarization tensor all neutrino- 
electron processes, allowed in a strong magnetic field, can be treated simultaneously. This calculation is accurate to 
first order in the Fermi theory, treats the effect of the magnetic field on the wave functions of the electrons exactly, 
and includes the effects of Pauli blocking. 

These calculations show that the rate of the standard unmagnetized processes may be modified greatly by the 
presence of a strong magnetic field, and that processes which only occur in a strong magnetic field are the dominant 
energy deposition processes for a range of field strengths and neutrino and antineutrino parameters. 

To make stronger conclusions about the physical implications of these processes in strong magnetic fields it is 
necessary to average the rates calculated here over thermal distributions of neutrinos appropriate to the astrophysical 
scenarios being considered. However, it is clear that there must be a reconsideration of the role strong magnetic fields 
and neutrinos may play in extremely energetic astrophysical events such as core collapse supernovae and gamma ray 
bursts. 
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